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■^ ■ ABSTRACT 

a'' 

I ■ We calculate the response of an ideal Michelson interferometer incorporating both dual 

recycling and squeezed light to gravitational waves. The photon counting noise has contri- 



^ ■ butions from the light which is sent in through the input ports as well as the vacuum modes 

at sideband frequencies generated by the gravitational waves. The minimum detectable 
gravity wave amplitude depends on the frequency of the wave as well as the squeezing 
and recycling parameters. Both squeezing and the broadband operation of dual recycling 
reduce the photon counting noise and hence the two techniques can be used together to 
make more accurate phase measurements. The variance of photon number is found to be 
time-dependent, oscillating at the gravity wave frequency but of much lower order than 
the constant part. 
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Laser-interferometric gravitational wave detectors [1] operate by sensing the difference 
in phase shifts imposed on the laser light in the two orthogonal arms of a Michelson type 
of interferometer by a gravitational wave. This phase shift manifests itself in the observed 
intensity change of the interference pattern. The sensitivity of the detector is determined 
by two fundamental sources of quantum mechanical noise : the photon counting error 
and the error originating due to fluctuations in radiation pressure on the mirrors. At the 
present level of the laser power available, the smallest detectable signal is limited by the 
photon counting statistics and various efforts have been made to increase the sensitivity 
level of these interferometers. 

Caves [2] first realized that the photon number fluctuations at the output could be 
understood due to the interference of the vacuum fluctuations of light which enters through 
the unused input port of the beam splitter with the ingoing laser light. He suggested that 
instead a squeezed photon state could be injected through this port to reduce the photon 
counting noise. For a Michelson interferometer operating on a dark fringe, most of the 
light escapes towards the laser source. Therefore, it had been suggested [3] that this light 
can be recycled by putting a mirror in front of the source to enhance the sensitivity of the 
interferometer. This technique is known as Power recycling. Brillet et al[4] argued that 
squeezing and power recycling are compatible with each other and that both can be used 
together to improve the signal-to-noise ratio. 

Gravitational waves modulate the phase of the laser light, thus generating sidebands 
which, travel towards the photodetector in an interferometer operating at the dark fringe 
[5]. These sidebands comprise the signal which can also be recycled by another mirror 
placed in front of the photodetector. The above method used in conjunction with power 
recycling is known as dual recycling [6]. 

It is, therefore, important to attempt an analysis of the quantum mechanical noise 
present in a dual recycling interferometer that also uses squeezed light and to investigate 
how well these two techniques work together. In this letter we report our results obtained 



for interferometers operating in the broadband mode. We arrive at a complete expression 
for the variance of the photon number fluctuations which is found to have a time-dependent 
component. The presence of sidebands significantly alters the noise. We calculate the 
minimum detectable gravitational wave amplitude as a function of its frequency as well as 
squeezing and recycling parameters and conclude that that the broad-band operation of 
dual recycling is compatible with the squeezed light technique and can therefore be used 
to enhance the sensitivity. 

We first evaluate the minimum detectable phase difference with both squeezing and 
dual recycling without considering gravitational waves. Referring to Fig.l, monochromatic 
light beams (of angular frequency ujq) in coherent and squeezed vacuum states enter ports 
1 and 2 respectively. The annihilation operators a and b represent light in the coherent 
and squeezed modes respectively. One may now write down equations for the intra-cavity 
electric field operators, Ea' and E^'. We assume that the distance between the recycling 
mirrors and the 50:50 beam-splitter has been adjusted such that Ea' and Eb' add in phase 
with the ingoing modes, Ea and Ei, respectively. We also assume that the beam-splitter 
introduces no phase shift upon refiection for a wave incident on the side of port 2 and a 
phase shift w for a wave refiected on the side of port 1. The quantities ti (t2) and ri (r2) 
represent the transmission and refiection coeff ! icients of the power (signal) rec 

One can obtain expressions for the annihilation operators a' and b' in terms of the 

input modes a and b. Then the annihilation operator describing the mode at the output 

of port 2 can be given as 

I 
Out2 :=t2b' -r2b= —[iatit2sme + b{4{cos6-ri) - r2M}], (1) 

where 

M = 1 + rir2 - (ri + ra) = (1 - ri)(l - ra) (2) 

and 9 is the phase difference of light between the two arms of the interferometer (at dark 
fringe, ^ = 0) . Then the mean and the rms value of the photon number at the output port 
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2 are found to be 
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respectively, where n is the mean number of photons in the coherent beam and r is the 
squeeze factor. In these expressions, we have neglected terms with coefficients sinh r since 
n ^ sinh r. At the dark fringe most of the laser light escapes towards port 1. However, 
for a very small phase shift 59, we obtain a very small change in the mean number of 
photons, 6N{6) at the output port 2. Equating the change to the rms value, we, therefore 
obtain the minimum detectable phase 56 at a dark fringe to be 
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As can be easily seen, for values of ri and r2 close to unity and a large squeeze factor, 
56 is considerably reduced. This shows that squeezing and recycling are compatible with 
each other and can be used together to increase the sensitivity of the interferometer. 

We now examine the case when a gravity wave of dimensionless amplitude, h{t) = 
ho sin Ugt, propagating along the z-axis impinges on an interferometer whose arms are 
oriented along the x and y axes. If the gravity wave interacts with the laser beam of 
frequency cuq, propagating along the y-axis for a time r then the phase picked up by light 
as a function of time is 

uoho 



mt) 



2uJn 



/ sm{ujgt) dt = eg sin ujg{t — —] 
Jt-T 2 



where uq is the laser light frequency and 
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(6) 



(7) 



Due to the quadrupolar nature of the gravity wave the phase acquired by the laser beam 
travelling along the x-axis is {—5(f)). The gravity wave thus modulates the phase of light 
in the two arms which gives rise to a time-dependent intensity [5]. 
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The positive frequency part of the electric field operator propagating along the y-axis 
can be written as 

where eo is the permitivity constant and Aq is the cross-sectional area of the quantization 
volume. One can now show [7] that after modulation, the positive frequency part of the 
electric field operator can be written as : 

„ , , f^ duj ( huJ \ • , . v^ / nUJa \l/2 ^ , . , 

E+{t)= iri^r^) ^ E {1 + ^)^ UeM^ + noOg), (9) 

Jo ^TT \2eoAoJ ^~^ '^o 



n= — cxD 



where a{uj + nujg) are the annihilation operators at newly-generated frequencies uj ± uJn 
and Jni^g) are the ordinary Bessel functions. 

So, for any wave originally present with frequency w, after modulation, one gets 
sidebands at a; ± nujg. Since Jni^g) ~ (^g)" fo^^ small €g and we are interested in terms 
of order 0{eg), we consider in our calculation only the first sideband on both sides. The 
sideband modes a{uj±ujg) and b{uj±ujg) which enter the interferometer through the ports 
1 and 2 respectively are all in their vacuum states. 

So, the intracavity electric field at the port 1 can be written as 



Ea'{t') 



{tiEa + riEa' + r2Ey + t2^b)e-'''^We-'^ 
+ (tiEa + riEa' - r2Eb' - t2Et)e'^'^^'^e'^ 



(10) 



where 9 is the constant phase offset between the two arms of the interferometer; t' = t+L/c 
and L is twice the arm length. The phase •0i is that acquired by light in traversing twice 
the path between the beam-splitter and the power recycling mirror. 

We can now write the Fourier transform of each electric field and then pull in the time 
dependent phase factor inside the integral to obtain 
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Now, we can easily see that the presence of e~*^'^*^*) factor inside each integral sign 
leads to phase modulation and subsequently to the generation of sidebands. We confine 
our attention only to the first sideband on both sides. The constant phase factor 9 = QAL 
is different for different frequencies, but, for simplicity, we assume it to be the same for all 
frequencies since the difference is very small {ug <^ ujq). We consider the coefficients ri, 
r2, ti, t2 to be independent of frequency. Since ujg <C ujq, we set the factor (1 + ujg/uj) -^ 1 
in Eq.(9). 

Now, one can write equations for the outgoing (primed) annihilation operators for 
Q = uj, uj ±ujg after dividing throughout by the same normalization constant. 



(12) 



where 



a'{n) exp[-z(OL/c + ^i(O))] =Ae-'^ + Be+'^ , 
b'{Q) exp[-z(OL/c + M^))] =Ae-'^ - Be+'^ , 

A =ti y^ a(0 + nujg)Jn{—h) + ri >^ a'(0 + nujg)Jn{—h) 

n n 

+ r2 2J&'(^ + nujg)Jn{-h) +t2 2J^(^ + nujg)Jn{-h), 

n n 

B =ti y^ a(0 + nujg)Jn{+h) + ri >^a'(0 + nujg)Jn{+h) 

n n 

- ^2 ^ h\Vt + nujg)Jn{+h) - t2 ^ 6(0 + nujg)J^{+h), 

n n 

where the index n can take values and ±1. The important quantity in the above equation 
is the phase factor (OL/c + 'i/'i(0)) that appears on the left hand side. The phase t/^i will be 
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(13) 



different for different frequencies O. However, the point to be noted here is that we adjust 
the distance between the recychng mirrors and beamsphtter in such a way that these phase 
factors become unity - a condition called 'on resonance'. This esssentially means that the 
laser light as well as the sidebands are resonant with the cavities formed by the recycling 
mirrors. This is termed as the broad-band operation of dual recycling. 

We now have six coupled equations for o! and V at three different frequencies (i.e. u) 
and uj ± ujg). One can arrange these equations in the following matrix form 

Pk^K = QkjAj, (14) 

where Pki and Qkj are two 6x6 matrices and 

A'- =Transpose(aQ, a'_, a^, ^q, 6'_, 6'_j_), 

(15) 
Aj =Transpose(ao, a_, a-|-,6o7 b-, b+), 

where (and from now onwards) indices 0, (— ) and (+) correspond to n = 0, —1 and +1 

respectively. So, for example, oq = a(a;o), ci- = a(i^o ~ ^g)^ ^+ = ^("^o + ^g) etc. 

So, all the six equations can be solved and the six primed annihilation operators 

can be written in terms of the six input (unprimed) annihilation operators through a 

6x6 matrix, P^i^Qki- If we define a^ and h^ as two 3x1 column vectors, i.e., a^ = 

Transpose(ao, a_, a+) and a^ = Transpose(ao, a_, a+), one can write the output fields at 

port 2 in a simple form 

=Xa^a^ + Yaubv, (16) 

where X^^. and Yau are two 3x3 matrices and (from now onwards) the greek indices take 
values 0, (— ) and (+). The values of different components of the 3x3 matrices Xq,^ and 
Yau are given below 

tit2 sin 6* 



-'^oo — X — — X_|__|_ — —i 



(l-ri)(l-r2)' 



^00 — Y — — Y+_|_ — 1, 



;i7) 



x+_ = x_+ = y+_ = Y. 



'+ 



+ --«+- 



-+ 



The components of order 0{eg) are 



0. 
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(18) 



''(l-ri)(l-r2) 

We essentially follow references [7,8] for the expression for the time-dependent pho- 
to cur rent, N . 



N 






X*X^pa}^ap + Y*Y^ph^J)p + X* Y^pa)j)p + Y^X^ph^^ap. 



(19) 



The mean number of photons, N is made up of a constant part, Nq = ^^ cj^Ca and the 
time-dependent part 5N{t) = '^^t^pc'l^Cis. The latter is essentially due to the beating of 
modes of two different frequencies which gives rise to the time-dependent part at cjg. There 
is also a time-dependent part at 2ujg but of order 0{e'^) and so we neglect it. Hence we get 

t?t| sin ^ 



^N>=no + SI{t) = N- 

(1 -ri)"=(l -r2)'^ 

The variance in photon number is given by 



[sin 6* + 4eg sincUgt]. 
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(21) 



(l-ri)2(l-r2)2 (l-ri)3(l-r2)3 

The terms appearing within the brackets of the constant part of Eq.(21) can be explained 
as follows: the first term is due to the coherent excitations (ao) superposed with coherent 



fluctuations as well as with the vacuum fluctuations of a_ and a+. The second term is 
coherent excitations superposed with vacuum fluctuations from 6_ and 6+, whereas the 
third term is due to the interference between squeezed and coherent light. The time- 
dependence of AN arises due to the beating of the time-dependent part of N with its 
constant part. Everywhere, we have neglected terms representing squeezed fluctuations 
being superposed on all the vacuum fluctuations since n ^ sinh r. The variance being 
time-dependent would mean that the frequencies separated by ujg are correlated although 
the spectrum is white. This has been referred to in the literature [9] as modulated shot 
noise. 

The minimum detectable gravity wave amplitude ho is now obtained by setting equal 
the maximum value of the 'signal', dN{t) to the maximum value of AN. This gives us a 
quadratic equation in ho 
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Since ho is already very small, we neglect the term ho/n. We flnally arrive at the expression 
for ho 



fto= '^' ,!,„.,+ (i^ii4lzi2)!+(i^ii)!(i^i£)!£:!:. (23) 

uoVfism{ujgT/2)]l 16 Stjtl tjtl 16 ^ ' 

At the dark fringe {9 = 0) the flrst term is negligibly small. For a large squeeze factor r and 
values of reflection coefficients ri and r2 close to unity, ho can be considerably reduced. 

If losses are introduced in the recycling mirrors, it would be possible to optimize ho 
in terms of the parameters r, ri, r2- Experimentalists usually implement internal [10], 
external [11] phase modulation (at MHz frequency) so that measurements of intensity can 
be made at sufficiently high frequency where the noise is really shot noise limited. Using 
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the equations described above, it should be possible to include both internal (or external) 
and gravity wave modulation (including losses) and calculate the minimum detectable 
gravity wave amplitude. Work in this direction is currently being pursued and will be 
communicated in future[12]. 

We are grateful to B.S. Sathyaprakash for his advice and many enlightening discus- 
sions. We also thank S.V. Dhurandhar and Kanti Jotania for suggestions. 
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Figure Caption 

Fig.l: A schematic diagram for dual recycling. BS- Beam-Splitter, EM- End Mirror, PRM- 
Power Recycling Mirror, SRM- Signal Recycling Mirror. 
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This figure "figl-l.png" is available in "png" format from: 
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